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ABSTRACT 

The ACS Nearby Galaxy Survey Treasury (ANGST) is a large Hubble Space Telescope (HST) 
Advanced Camera for Surveys (ACS) treasury program to obtain resolved stellar photometry for 
a volume-limited sample of galaxies out to 4 Mpc. As part of this program, we have obtained deep 
ACS imaging of a field in the outer disk of the large spiral galaxy M81. The field contains the outskirts 
of a spiral arm as well as an area containing no current star formation. Our imaging results in a color- 
magnitude diagram (CMD) reaching to iriF8i4iy = 28.8 and f&fw6W = 29.5, one magnitude fainter 
than the red clump. Through detailed modeling of the full CMD, we quantify the age and metallicity 
distribution of the stellar populations contained in the field. The mean metallicity in the field is 
— 1 < [M/H] < and only a small fraction of stars have ages £j 1 Gyr. The results show that most 
of the stars in this outer disk field were formed by z ~ 1 and that the arm structure at this radius 
has a lifetime of ^ 100 Myr. We discuss the measured evolution of the M81 disk in the context of 
surveys of high-redshift disk galaxies and deep stellar photometry of other nearby galaxies. All of 
these indicate that massive spiral disks are mostly formed by z~l and that they have experienced 
rapid metal enrichment. 

Subject headings: galaxies: individual (M81) — galaxies: stellar populations — galaxies: spiral — 
galaxies: evolution 
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1. INTRODUCTION 

Analytic and numerical models indicate that spiral 
disks should grow and evolve with time due to con- 
tinued gas accretion, interactio ns, spiral density waves , 
and internal vis c ous e v olution (|Fall fc Efstathioul 1 19801 
Dalcanton ct al 
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straints on the evolution of disks have come largely 
through identifying changes in the bulk properties of the 
galaxy population from in situ measurements at high red- 
shifts (up to z ~ 1 — 1.5). However, these observational 
attempts t o confirm disk evolution models give conflict- 
ing results (Simard ct al. 1999; Ravin dranath et all2004t 
Truiillo et al.ll20 04: Bar den et al. 2005: M elbourne et all 
2007t ICameron fc Driver! I2007D . likely due to selection 
effects that are difficult to quantify. Fortunately, the 
evolution of disks can be independently constrained with 
photometry of the resolved stellar populations in nearby 
galaxies. Such photometry provides a fossil record of the 
formation and evolution of the disk and complements the 
findings of redshift surveys. 

Resolved stellar photometry provides the most detailed 
data with which to determine the star formation his- 
tory (SFH) of a galaxy. By fitting stellar evolution 
models to an observed color-magnitude diagram (CMD), 
we can recover the stellar ages and metallicities that 
best reproduce the color and magnitude distribution of 
a galaxy's stars. To fully tap this capability, the ACS 
Nearby Galaxy Survey Treasury (ANGST) has under- 
taken a program to measure resolved st ellar photometry 
for a volume- limited sample of galaxies ([Dalcanto n et all 
2008). Within this volume, large galaxies dominate the 
stellar mass. Of these, the largest and most well-studied 
at all wavelengths is M8 1, a massive SA(s)ab s piral disk 
at a distance of 3.9 Mpc (jTikhonov et al.ll2005f) with low 
foreground extinction (Av =0.27: ISchlegel et al.lll998h . 



-..20061) places 
23: iKochanek et al.ll200l[ ). 



Its luminosity (M#-=-24; Skrutski e et al 
it at 2.5L* (assuming Mk 
making it characteristic of the disk galaxies seen in red- 
shift surveys out to z ~ 1 (e.g.. lOvaizu et al1l2008t ). 
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M81 provides a key laboratory for using resolved stel- 
lar photometry to look for relics of the evolution seen 
at high-redshift. Several surveys suggest that large disk 
galaxies like M81 have their disks in place by z ~ 1, 
after which luminosity evolution dominates (|Lilly et ID 
19981: iRavindranath et all 12004 iPapovich et all 120051 : 
Sargent et al l 120071 : Mel bourne et al.ll2007l ). In this con- 
text, understanding the SFH of nearby large disk galax- 
ies like M81 provides complementary insight into their 
evolution that is free of any biases contained in redshift 
survey samples. 

Much work has been done to unde rstand the evo- 
lution of M81 from its star clusters (|Ma et all 120051 : 
iPerelmuter et all 1 1995ft . its X -ray source population 
(|Swartz et all 12003ft . and its young supe rgiant stars as 
seen in the near infrared (Davidgc 2006). Its structure 
and evol ution have been studied in detail from ultraviolet 
to radio (iPerez-Gonzalez et alj|2006t IWillner et al.ll2004t 
IGordon et ail f2004: iLi et al. 200J" lWestpfahlll998h . but 
very little work has been done to integrate these results 
with the resolved stellar populations. 

Most work on resolved stellar populations in M81 has 
relied on the Hubble Space Telescope (HST) to resolve the 
individual stars, as only stars brighter than the red gi- 
ant branch (RGB) t i p are resolved from the ground (e.g., 
iMadore et ail 119931) . iTikhonov et all (|2005| ) studied the 
resolved stellar populations of M81 with archival WFPC2 
and ACS data, and there are several HST programs 
that are curr ently underway with greater depth and spa- 
tial coverage (Jde Jong et aLll2007tlSaraied ini 2005; Zezas 
120051: |Huchra|l200|). The results from these programs 
will allow detailed comparisons with those from inte- 
grated light studies. 

Herein we describe the ANGST measurements of the 
SFH of the outer disk of M81 using F606W and F814W 
stellar photometry. The outer disk has low enough 
crowding that precision photometry can be obtained be- 
low the red clump with ACS, and provides leverage with 
which to potentially constrain size evolution of the galac- 
tic disk. §2 details our data acquisition and reduction 
techniques. §3 discusses our analysis of specific portions 
of the color-magnitude diagram (CMD), and §4 details 
our analysis of the full CMD to determine the SFH of the 
field. Finally, §5 summarizes the conclusions drawn from 
these results, placing our measurements in the context 
of disk galaxy surveys and resolved stellar populations of 
other large galaxies. We adopt a cosmology with H = 73 
km s" 1 Mpc-\ fl A = 0.76, and tt M = 0.24 for all look- 
back time calculations, and we assume a distance to M81 
of 3.9 Mpc for conversions of angular measurements to 
physical distances. 

2. DATA ACQUISITION AND REDUCTION 

From 2006-Nov-16 to 2006-Nov-22, we observed a field 
in the outskirts of the M81 disk located at R.A. (2000) 
= 148.644625 (09:54:34.7), decl. (2000) = +69.2804 
(+69:16:49.4) with a rotation angle PA_V3=89.81. Fig- 
ure Q] shows an outline of the field location, which is 14' 
(16 kpc at M81) out along the major axis and corre- 
sponds to 5 scale lengths an 18 effectiv e radii of the 
bulge (h r ~ 3 kpc and bulge R e ~ 0.9 kpc, ' Kendall et al.l 
2008). The equivalent location in M31, as shown on the 
inset in Figure [TJ suggests that the disk population is 
likely to still dominate at this radius if M81 and M31 



are similar. T his suggestion i s cons istent with the ongo- 
ing analysis of lde Jong et al.l (|2007| ). which finds that the 
halo population does not dominate until galactocentric 
distances ^ 20 kpc. 

We obtained 9 full-orbit exposures with the ACS 
(|Ford et all [1998ft through the F606W (wide V) filter, 
and 11 full-orbit exposures through the F814W (I equiv- 
alent) filter. These data totaled 24132 s and 29853 s 
of exposure time in F606W and F814W, respectively. 
Routine image calibration, including bias corrections and 
flat-fielding, were performed by the HS T pipeline, OPUS 
version 2006_5, CALACS version 4.6.1. We processed 
the imag es by running the mul tidrizzle task within 
PyRAF (|Koekemoer et al.l 12002). This procedure was 
used only to flag the cosmic rays in the individual im- 
ages, after which, photometry was measured simultane- 
ously for all of the objects in the unco mbined images 
using the software package DOLPHOT (iDolphinl [2000h 
including the ACS module. This package is optimized 
for measuring photometry of stars on ACS images using 
the well-characterized and stable point spread function 
(PSF), calculated with TinyTim. 14 The software fits the 
PSF to all of the stars in each individual frame to find 
PSF magnitudes. It then determines and applies the 
aperture correction for each image using the most iso- 
lated stars, corrects for the charge transfer efficiency of 
the ACS detector, combines the results from the individ- 
ual exposures, and converts the measured count rates to 
the VEGAmag system. 

The DOLPHOT output was then filtered to only al- 
low objects classified as stars with signal-to-noise >6 in 
both filters. The list was further culled using sharp- 
ness (\F606W sharp + F8UW sharp \ < 0.27) and crowd- 
ing {F606W crowd + F8UW cr0 wd < 0.1). The sharpness 
cut was chosen based on the distribution of values in the 
original catalog. This distribution is shown in Figure [2] 
and flattens at a value of ~0.27. The crowding parame- 
ter gives the difference between the magnitude of a star 
measured before and after subtracting the neighboring 
stars in the image. When this value is large, it suggests 
that the star's photometry was significantly affected by 
crowding, and we therefore exclude it from our catalog. 
Quality cuts based on the x values were also considered, 
but they were rejected when a correlation was found be- 
tween x and the local background. Our final star cata- 
log for the field contained 120912 stars detected in both 
F606W and F814W, and the resulting CMD is shown in 
Figure [H 

The same software package was used to perform arti- 
ficial star tests using identical measurement techniques 
and quality cuts. We ran 2.5xl0 6 artificial stars to char- 
acterize our photometry errors and completeness as a 
function of color, magnitude, and position. In each it- 
eration, a single star was added to the images, and the 
photometry of the images was remeasured in the area 
where the star was added, including a radius of the PSF 
size plus the photometry aperture size to include the pho- 
tometry of all stars whose photometry could be affected 
by the existence of the artificial star. Half of the artifi- 
cial stars were sampled randomly in color and magnitude 
covering the full range present in our observed CMD plus 
an additional magnitude fainter to account for upscatter 

14 http:/ /www. stsci.edu/software/tinytim/ 
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of faint stars into our recovered magnitude range. The 
other half were sampled following the color and magni- 
tude distribution of our observed CMD after extrapo- 
lating the distribution to fainter magnitudes to account 
for upscatter. The artificial stars were distributed ran- 
domly over the field of view. The photometric errors 
measured from our tests are shown in Figure [51 and 
the completeness measured from our tests is shown in 
Figure O When fitting the color-magnitude distribu- 
tion of the stars, we included only stars brighter than 
a 60% completeness limit as measured from the artifi- 
cial star tests (m F eoew = 29.1 and m F si4w = 28.4). At 
this depth and this Galactic latitude (b = 40.9°), the 
expected number of Galactic foreground stars is ,$ 20 
axcmin , suggesting foreground contamination in our 
field of ;$ 200 objects or <0.2% of the total number of 
stars. 

3. RED CLUMP AND ASYMPTOTIC GIANT BRANCH 
BUMP ANALYSIS 

3.1. Overview 

Before attempting to understand the complexities of 
the CMD as a whole, we first focus on a few key stellar 
evolution features in our full-field CMD that help to give 
a broad sense of the range of ages and metallicities of 
the stars in the field (Figure [3]). The vertical plume at 
mF606W — 1TIF814W = is due to the young population 
of stars that is still on the main sequence. The handful 
of stars populating the diagram brightward and slightly 
redward of this plume are massive core He-burning stars, 
which are young stars that have recently evolved off of 
the main sequence. We plot the spatial distribution of 
the main sequence stars in Figure [4j the young stars are 
clearly concentrated in the inner disk (seen in the bottom 
of the image) and in the extension of the spiral arm seen 
in Figure [1] 

Several older populations are seen as well. The dense 
clump of stars at mpsi4w ~ 27 . 8 is the red clump 
(jCannonl [l970t lAlves fe Saraiedinil 1 19991) . which corre- 
sponds to the core-He burning phase of all intermediate- 
age and old populations which are neither too old, 
nor too metal-poor, to develop the horizontal branch. 
Brightward of this feature is another, less prominent 
peak, known as th e asymptotic giant branch (AGB) 
bump (|Gallartlll998| ). which corresponds to the so-called 
early-AGB phase of low-mass stars, during which He 
shell burning transits from a very extended to a thin 
shell. Extending vertically through both the red clump 
and AGB bump is a broad RGB with a well-defined tip, 
indicating the presence of a range of ages and metallici- 
ties in the field. Brightward of this tip is a relatively small 
number of thermally-pulsating AGB stars, i.e. stars in 
the stage of double shell burning which undergo recur- 
rent He shell flashes. The most luminous among these 
thermally-pulsating AGB stars have intense winds that 
eventually shed their outer layers to become planetary 
nebulae. 

3.2. Red Clump and AGB Bump Fitting Method 

To investigate the general characteristics of the stel- 
lar populations in our field, we first performed precision 
measurements of the very well-defined red clump and 
AGB bump, seen at mpsuw ~ 27.8 and iriFsuw ~ 26.7 
in the diagrams in Figure [3] respectively. These are areas 



of the CMD containing a high density of stars, allow- 
ing their precise locations to be measured and directly 
compared to stellar evolution models. The colors and lu- 
minosities of these features have be en used to constrain 
the properties of ste llar populations (Rei kuba et al.ll2~005t 
iTanaka et aTl I2008T ). We include such a measurement 
here to take advantage of the depth and quality of our 
photometry and to allow for intercomparisons with other 
such analyses. However, as we discuss below, uncer- 
tainties with the models of these discrete features sub- 
stantially limit the interpretation of these measurements, 
making simultaneous fitting of the full CMD preferable; 
we present such fits in § 4. 

We measured the magnitudes of the red clump and 
AGB bump by creating a magnitude histogram from our 
photometry and fitting it with a combin ation of a line 
and 2 Gaussians, following the methods of lReikuba et al.1 
(2005). We then fit Gaussians to histograms of slices 
through color space taken at the measured peak mag- 
nitudes (Figure [7|). The red clump in this outer re- 
gion of the M81 disk has m FSUW = 27.792 ± 0.002 
(FWHM=0.450) and m F606W - m F saw = 0.855 ±0.001 
(FWHM=0.363). The AGB bump has m FSUW = 
26.72 ± 0.01 (FWHM=0.27) and m F60 6W ~ m F8 uw = 
0.974 ± 0.002 (FWHM=0.346). These are the raw mea- 
surements for these features; no extinction correction has 
been applied. 

We then measured the locations of model red clumps 
and AGB b u mps by populating the isochrones of 
iGirardi et all (120021) inclu ding updates of AGB mod- 
els in iMarigo et al.l (12008ft . 15 We used the StarFISH 
(|Harris fc Zaritsky|2001f ) task testpop, to produce a grid 
of CMDs for discrete ages and metallicities, assigning ex- 
tinctions, distances, errors, and completeness measured 
from our data (see §S]). We then ran our red clump and 
AGB bump fitting technique on this grid of model CMDs 
to measure the location of these features as a function of 
age and metallicity. 

We determined the model red clump and AGB bump 
with F606W and F814W magnitudes closest to those 
observed in our field, assuming an adopted distance of 
(m - M) = 27.93 and extinction of A v = 0.27. These 
distance and extinction values resulted in a grid of model 
red clump and AGB bump values that did not intersect 
with our measured values. We therefore show a model 
grid with a farther assumed distance ((m — M)q = 28.0) 
and higher extinction (Ay — 0.3) in Figure El to pro- 
duce a grid that does intersect with our observed val- 
ues. From this analysis, the models suggest that the red 
clump is dominated by stars with an age of ~2-3 Gyr 
and [M/H]~-0.4 and that the AGB bump is dominated 
by older, more metal-poor stars with an age of ~10 Gyr 
and [M/H]— 0.7. 

Because the precise colors and magnitudes of the red 
clump and AGB bump do not provide reliable conclu- 
sions (see S I3.3|) . we also investigated the relative num- 
bers of stars contained within the red clump and AGB 
bump features. This measurement should be insensitive 
to distance and extinction. We calculated the integral of 
the functional fit to the magnitude histogram including 
and excluding the Gaussian terms corresponding to each 
feature. The difference between these integrals provided 

15 http://stev.oapd.inaf.it/cmd 
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the number of stars within each feature. The ratio of 
the number of stars in the AGB bump to that contained 
in the red clump {N AGBb : N RC ) was 0.044±0.008. Er- 
rors are the standard deviation of the same measurement 
made on 100 random samples of our data and are dom- 
inated by the small number of stars in the AGB bump. 
Comparing this range of ratios to the ratios obtained by 
running the same calculations on the functional fits to 
models suggests that the dominant stellar population of 
our sample is metal-rich and older than 3 Gyr (see Fig- 
ure [9]) . The ratio is therefore consistent with the color 
and magnitude of the model AGB, but less so with those 
of the red clump, in that the model red clump is fainter 
than it should be for a population old enough to have 
this fraction of AGB bump stars. 

3.3. Interpretation of Red Clump and AGB Bump 
Fitting Results 

Although analyses similar to those above are becoming 
common, the red clump and AGB-bump colors, magni- 
tudes, and relative numbers cannot provide conclusive 
results about the age and metallicity distribution of the 
population. The reasons for the inconclusive results in- 
clude the sensitivity of the results to known deficiencies 
of stellar evolution models as well as to the assumed dis- 
tance and reddening of the stars in the CMD. 

We found our fits to the color and magnitude of the 
red clump and AGB bump difficult to interpret because 
they are particularly sensitive to weaknesses of the stel- 
lar evolution models. The majority of the red clump 
models are faint compared to our observed red clump. 
While the trends of the colors and luminosities of the red 
clump and AGB bump with age and metallicity (older 
and more metal-rich are fainter and redder) are robust 
against model updates and different model prescriptions, 
the absolute colors and luminosities are less stable. For 
example, red clump model luminosities depend primar- 
ily on opacities and neutrino losses during the previous 
RGB phase, which ar e still subject to sig nificant varia- 
tions in the literature. iGirardi et alj (|2002f ) models of the 
red clump are known to be some of the faintest because 
of th eir treatment of these processes (fCastcllan i et al.1 
2000). This effect can be seen in Figure [8j where nearly 
all of the model red clumps are fainter than our observed 
red clump even though most of the model AGB bumps 
are brighter than our observed AGB bump. 

The NAGBb ■ Nrc ratio is also difficult to interpret. 
The ratio predicted by the models is sensitive to the as- 
sumed helium content and the treatment of convective 
cores, including complex processes such as overshooting 
and semiconvection, rendering any conclusions based on 
this ratio alone unreliable. 

Furthermore, fits to the color and magnitude of the 
red clump and AGB bump depend on the extinction 
and distance values applied to the models, as changes 
of only —0.1 mag in color or luminosity of the fea- 
tures correspond to differences of ~0.3 dex in metal- 
licity and log(age). The effects of distance and red- 
dening uncertainty are shown with the arrows in Fig- 
ure [5] While the adopted distance affects mainly the 
measured age, the adopted reddening largely affects the 
measured metallicity. More specifically, if we assume 
only foreground extinction value (Ay — 0.27), our best- 
fitting ages remain the same, but the metallicities become 



[M/H]~-0.4 and [M/H]~- 0.1 for the AGB bump and red 
clump, respectively. If we assume a distance modulus of 
(m — M)o = 27.9, the best-fitting ages and metallicities 
are [M/H]~- 0.7, [M/H]~- 0.5 and -13 Gyr, -5 Gyr for 
the AGB bump and red clump respectively. 

In summary, while studying the red clump and AGB 
bump properties in detail is helpful for getting some sense 
of the overall age and metallicity of the population, or 
perhaps for constraining relative ages and metallicities 
among different galaxies, the results' sensitivity to uncer- 
tainties in stellar models, distance, and extinction, make 
this analysis less than optimal for obtaining reliable in- 
formation about the stellar populations. Furthermore, 
no single age and metallicity will precisely fit the ob- 
served values because these features contain populations 
of a range of ages and metallicities. These difficulties 
with single component fitting make it necessary to per- 
form more sophisticated statistical fits to the entire CMD 
to decipher the range of ages and metallicities present in 
our field. Full CMD fitting helps to reduce the effects 
of model deficiencies, providing an overall picture of the 
age and metallicity of the population even if some details 
of the models are wrong. 

4. FULL CMD FITTING 

4.1. Fitting Technique and Models Used 

We measured the complete SFH, including the star for- 
mation rate and m etallicity as fun ctions of age, using the 
MATCH package (jDoIphinl 12002] ). This software fits the 
entire observe d CMD by populati ng the stellar evolu- 
tion models of IGirardi et al.l (I2002L with updated AGB 
models from lMarigo et al.ll2008f ) with a given initial mass 
function (IMF), finding the distance modulus, extinction, 
and linear combination of ages and metallicities that best 
fit the observed color and magnitude distribution (see de- 
tails in iDolphinl l2002f) . 

As was the case for the simplified AGB bump and red 
clump analysis, the results of full CMD fitting rely on 
the stellar evolution models used to fit the CMD. There 
are known discrepancies between different sets of stellar 
evolution models in the RGB, red cl ump, and AGB bum p 
phases of evolution (see Figure 10 of lGallart et aD l2005l. 
The reasons for the discrepancies are many, from the 
adopted input physics to the uncertainty in mass loss 
during the RGB and AGB phases of evolution. These 
differences affect the location and morphology of these 
regions of the observed CMD. 

On the other hand, the shape of the red clump re- 
gion does contain information about the SFH. The fit- 
ting algorithms are sensitive to this information if they 
are set to fit the distance and reddening values indepen- 
dently in order to compensate for offsets between model 
red clump positions. For instance, the RGB colors are 
sensitive to both age and metallicity, whereas the RGB- 
red clump color diffe rence is more sensitive to the age 
(|Hatzidimitrioul[l991h . While many age- metallicity mod- 
els could fit the overall RGB color distribution, when 
the red clump is present its color difference forces the 
SFH-recovery towards the best-fitting age distribution, 
provided that it is allowed to fit also the red clump mag- 
nitude (e.g. by adjusting the distance). The final result 
is that this area of the CMD provides a reliable SFH, 
even if the models for the red clump contain offsets. To 
test this point, we ran a different fitting code and stel- 
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lar evolution library to fit our full-field data. We ap- 
plied t he IAC-STAR/IAC-POP/ MINNIAC CMD fitting 
codes (jAparicio fc Hidalgdl200l 120081 Hidal go, S. L. et 
al. 2008, AJ, in preparati on) to our data using the BaSTI 
stellar evolution models (|Pietrinferni et all 12004). The 
broad trends on the resulting SFH were totally consis- 
tent. 

Thi s point is strengthened by the tests of lBarker et al.l 
(|2007l ). who compared results for SFHs as determined 
with different sets of stellar evolution models and found 
differences in the details but agreement for general 
population characte ristics . Furthermore, t he results of 
iTolstov et all (|f998l ) and iCole et all (|2007t ) show that, 
while SFHs based on the main-sequence turnoff pro- 
vide more reliable age information at older ages than 
those from shallower photometry, the general trends de- 
termined in shallower photometry are robust. To avoid 
drawing conclusions based on the finest details in age and 
metallicity applied to fit the CMD (see §[13]), we bin our 
results to an age and metallicity resolution where both 
the CMD fitting method and stellar evolution models are 
well-tested and the limitations are known. 

4.2. Field Division 

Our field contained a portion of a spiral arm running 
through the north and east portion of the image. The 
bifurcated arm, which may have been created by an inter- 
action, can be seen in the locations of the main-sequence 
stars, which are mainly limited to the northeast portion 
and the southeast corner (Figure 2]) . This distribution 
shows that the northeast spiral arm is split, with a spur 
to the nor th, consistent with H I m aps of lYun et all 
(|1994h and lAdler fc Westpfahll (|1996t ). The recent star 
formation in this structure is likely to be at least partially 
due to a tidal stream from one of several rece nt interac- 
tions with other galaxies in the M8f group (|Yun et al.l 
11994 lYur]|f999D . The southeast region of the image is 
more crowded as well, where it skims the inner disk. 

For measuring the SFH, we divided the field into re- 
gions inside the arm, outside the arm, and the more 
crowded inner disk, as shown in Figure |4j Since the 
crowding and extinction are different in each of these 
subregions, we determined the error and completeness 
characteristics separately for each subregion. The final 
CMDs for the arm and interarm subregions are shown 
along with the full CMD in Figure [3] 

4.3. Fitting Parameters 

To model a full CMD, several fitting parameters must 
be chosen. These choices include the binary fraction, 
IMF slope, the area of the CMD to include in the fit, the 
approximate distance and mean extinction to the stars in 
the field, and the binning of the stellar evolution models 
in time and metallicity. Below we discuss how we chose 
these parameters and how the choices impact our results. 

When populating the mo del isochrones, w e assumed a 
binary fraction of 0.35 and a lSalpeterl (|1955l) IMF. As has 
been shown by other studies using this technique (e.g., 
I Williams et"aLll2007t iBarker et al.ll2007f h for photometry 
that does not include the main sequence of old stars, the 
IMF assumed does not affect the relative star formation 
rates in the SFH, but does affect the normalization of the 
SFH. This normalization effect occurs because the red 
clump, RGB, and AGB probe a narrow range of initial 



mass. Therefore changing the slope has little effect on 
the relative star formation rates, but a large effect on the 
extrapolation of the mass contained in the underlying 
unresolved low-mass main-sequence stars. Since we are 
not attempting to determine the precise star formation 
rate but rather are interested in the relative SFH within 
the field, a Salpeter IMF is sufficient for our purposes. 

When fitting the color-magnitude distribution of the 
stars, we included all stars brighter than our 60% com- 
pleteness limit for the full field to avoid large complete- 
ness corrections in our model fitting. This completeness 
cut corresponds to mF606W — 29.1 and ttifsiw = 28.4, 
and only stars brighter than these magnitudes were in- 
cluded in our model fits. We allowed the distance mod- 
ulus to range from 27.75 to 28.2 and the extinction to 
range from A v = 0.1 to A v = 0.7, which allows MATCH 
to determine the systematic errors that result from small 
changes to these parameters and to optimize the overall 
CMD fit even in the presence of localized deficiencies 
in the model isochrones. While differential extinction 
can, in principle, be a problem for fitting a CMD with 
a single extinction value, our M81 field lies in a region 
with very little visi ble dust structure in Spitzer maps 
(|Kendall et al.ll2008| ). Furthermore, modest amounts of 
differential extinction have been shown to have little im - 
pact on the results of full CMD fitting (| Williams! 120021 ). 

The best- fitting distances (see Table [1| were all consis- 
tent with (to — M)o — 27.9 within the errors, which agree 
with the Cepheid dista nce from the HST key project 
((to - M) = 27.8 ± 0.2; iFreedman et al.lll994[ ) and the 
distance determined by measuring the tip of the RGB in 
archival WFPC2 pho tometry ((to - M) = 27.93 ±0.04; 
iTikhonov et~aTll2005h . 

The best-fitting extinction for the interarm region was 
consistent with the value of Ay=0.27 obtained from 
ISchlegel et al.l (|1998f ) for the Milky Way foreground ex- 
tinction (see Table [T|) . We note that the arm subregion 
had a measured extinction value that was significantly 
higher than the Milky Way foreground, reflecting the 
higher dust content expected in the arm region. 

We used a fine logarithmic time and metallicity reso- 
lution (0.1 dex) when fitting the CMD to allow the best 
possible fit to the data. The fit to the full data set is 
shown in Figures [TO] and QT] The full CMD fitting con- 
firms that the models are not able to perfectly reproduce 
the red clump and AGB bump, as expected from our dis- 
cussion of uncertainties in § [3] These are the only two 
features of the CMD that show significant differences be- 
tween the data and the best-fitting model (bottom-right 
panel of Figure [TO]) , confirming that the models of these 
features still need improvement. After performing the 
full CMD fit, we binned the age and metallicity results to 
coarser time resolution to reduce our SFH errors. There- 
fore while our fit did not force the star formation rate or 
metallicity to be constant within a given temporal bin, 
our final SFHs only show the mean rate and metallicity 
within each temporal bin to avoid drawing conclusions 
based on details of the fit that may not be robust against 
changes in models and fitting methods. 

4.4. Error Determination 

In addition to measuring the most likely SFH for our 
field, we ran Monte Carlo tests to determine the ran- 
dom uncertainties of the fits. We generated CMDs by 
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randomly drawing stars from our observed CMD, allow- 
ing each star to be drawn any number of times. We 
then measured SFHs for the resulting CMDs and calcu- 
lated differences from the best-fit to the actual data. We 
generated f 00 samples with the same number of stars as 
our observed CMD, adopting each of the SFHs measured 
from our data to determine the random errors for each 
subfield. These errors were then added in quadrature to 
the systematic errors, determined by fitting the CMDs 
with a range of possible distance and reddening values, 
to provide the our final errors on the rate, metallicity, 
and cumulative fraction of stars formed as a function of 
time. 

We note that our choice of time bin widths was 
very coarse for ages >2 Gyr to reflect our sensitiv- 
ity to age from the RGB and AGB. From ~2-8 Gyr 
(9.3<log(age)<9.9), the age distribution comes mostly 
from the relatively small number of stars on the AGB, 
which is a very short-lived and difficult to model stage 
of evolution. From 8-14 Gyr (9.9<log(age)<10.15) age 
has very little effect on the morphology of the AGB and 
RGB features. In addition, our metallicity measurements 
at young ages (<100 Myr) have large errors. Our only 
metallicity information for these ages comes from the 
relatively small number of stars on the short-lived He- 
burning sequences. Despite these unavoidable sources of 
uncertainty in age and metallicity, overall we are able 
to obtain very reliable estimates of the relative contribu- 
tions of stars of old (>8 Gyr), intermediate (2-8 Gyr), 
and young (<2 Gyr) ages. Furthermore, we obtain re- 
liable metallicities covering all but the youngest ages. 
Finally, the cumulative age distribution is stable against 
uncertainties at intermediate ages because systematic er- 
rors in the star formation rates in adjacent time bins are 
typically anti-correlated. 

4.5. The SFH of the Full Field 

The SFH of the entire field and the three independent 
subfields are shown in Figures [T2T[T71 For the history of 
star formation in the full field (Figure [12] and Table [2]) , 
we find that more than 50% of the stars currently at 5 
scale lengths from the galactic center formed by z ~ 1, 
and ~70% formed by z ~ 0.5. The bulk of the stars have 
— 1 [M/H] ^ 0, with no significant metallicity differ- 
ence between the arm and interarm regions and with no 
evidence for a significant metal-poor component. Such a 
me tallicity is consis t ent w ith the mean metallicities found 
bv iTikhonov et alJ (|2005l <[M/H]> 0.65) from stel- 
lar photometry in 5 WFPC2 fields spanning a range of 
galactocentric distances out to ^5 scale lengths. It there- 
fore appears that M81 was chemically enriched very early 
in its history. Moreover, the metallicity has not changed 
by more than ~1.0 dex (and possibly as little as 0.5 dex) 
since its very early history. 

Although our metallicity results for recent times have 
large errors, their values (-0.5< [M/H] <0.0 at ages from 
10 Myr to 100 Myr; Figure [P2]) are consistent with the 
gas-phase metallicity at this radius ([O/ H] ~ —0.3) in- 
ferred from the abundance gradients of IZaritskv et al.l 
(1994), suggesting that the gas responsible for this struc- 
ture is from M81 and not recently accreted gas from 
a more metal-poor interacting satellite. However, this 
metallicity does not rule out the possibility that the 
gas came from a metal-rich satellite, such as NGC 3077 



([M/H]~0, lMartinlfT99l . 

4.6. The SFHs of the Subregions 

Figure [TBI shows the SFHs of the three subregions over- 
plotted with the SFH of the total field. It is clear that 
the sum of the results of the three subregions measured 
independently is consistent with the result of the entire 
field measured at once. The consistency between the 3 
independently measured SFHs and the total SFH for the 
field confirms that our measurement technique provides 
reliable and consistent results in a field containing re- 
gions with moderately different crowding and extinction 
properties. Furthermore, the more crowded portion of 
the field, closest to the galaxy center, has a SFH more 
similar to that of the interarm region than to that of the 
arm region, showing that this area is not part of an arm 
and is of higher density only because it is closer to the 
galaxy center. 

While these three regions should have obvious differ- 
ences in their recent SFHs (based on the distribution of 
main sequence stars; cf. Figured]), the older stellar pop- 
ulations should be well mixed on timescales longer than 
the dynamical time ( £ 0.5 Gyr). Indeed, we measure 
similar SFHs at times 0.5 Gyr. At younger ages, the 
interarm region shows a clear deficit of stars younger than 
a few hundred Myr. Fractionally, this difference corre- 
sponds to ^2% of the total star formation (see Figure fl"7l 
and Table [2]). The arm region of our field contains es- 
sentially all of the stars with ages ^ 100 Myr. This 
shows that the stars have not scattered out of the arm 
structure on this timescale. In addition, there is a sig- 
nificantly higher fraction of stars with ages ,$,0.5 Gyr 
in the arm region than in the other regions. It there- 
fore appears that the arm structure in this region has 
a lifetime of £ 100 Myr and that full dispersion occurs 
on a timescale of ~0.5 Gyr. The characteristic width of 
the arm region, divided by the typical lifetime gives an 
approximate speed for the diffusion of star s. Assuming 
an arm width of -1 kpc (|Westpfahll fl99l divided by 
a lifetime of ~100 Myr yields a diffusion speed of ~10 
km s _1 , which is similar to th e velocity dispersion of B 
stars in the Milky Way disk (jDehnen fc Binnevl[l998l ). 
Therefore the timescales we are measuring for this arm 
feature are consistent with the stellar kinematics of our 
own Galactic disk. Since our data contain main-sequence 
and He-burning sequence information for the populations 
of these ages, spatially resolved SFHs, with a time res- 
olution of ~ 25 Myr over the last ~ 300 Myr are cur- 
rently in progress and will be presented in future papers 
(|Gogarten et all"2008t B. Williams et al., in preparation). 

5. CONCLUSIONS 

We have performed deep resolved stellar photome- 
try of an ACS field in the outer disk of M81. The 
metallicities of the stars in the field appear to have in- 
creased by at most ^0.5 dex over the past 10 Gyr from - 
1.0 < [M/H] < -0.5 to -0.5 < [M/H] < 0, suggesting early 
enrichment and a mechanism at work diluting the enrich- 
ment products of roughly continuous star formation. 

Similar behavior has been seen in the thick disk of 
the Milky Way and several other nearby large galaxies. 
Figure [18] shows the metallicities of the dominant stel- 
lar population for several nearby galaxies and the Milky 
Way thick disk as a function of galaxy luminosity, field 
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location, and dominant stellar age. The Milky Way thick 
disk has characteri stic ages of 9-12 Gyr and metallicities 
of fM /Hl > - 0.7 (tGlimore et al J 11993 iProchaska et all 
120001 lAllende Prieto et al.ll2006fh Thick disks in some- 
what smaller galaxies also appear to have metallicities 
>-1.0 (jSeth et al.ll2005f ). In addition, deep iJSTphotom- 
etry of the M31 disk has revealed populat ions dominated 
by stars in this ag e and metallicity ra nge (Willi Mn^[200l 
lOlsen et alj|2006t iBrown et al l 120061), as has deep pho - 
tometry of the outer disk of M33 (|Barker et all 120071) . 
All of these results suggest that the histories of these 
large disk galaxies may have been similar. 

Furthermore, even in the outskirts of large ellipti- 
cals, such as NGC 5128 ([M/H] ~-0 6 and age~8 Gyr 
iReikuba et all 12001 and NGC 3377 (|Harris et al.ll2007L 
[M/H] ~-0.6 and age >3 Gyr,), the age and metallicity 
of the dominating population is similar to that seen in 
the outer regions of M81 and other large disks. Taken 
together, these data point to a large galaxy formation 
scenario with rapid early enrichment to [M/H] ^-1.0, 
before lookback times of ~7 Gyr (z > 0.9). 

There is also significant evidence that most of the 
stars in this field were formed by z ~ 1. Such a re- 
sult is in general agreement with the findings of sev- 
eral recent galaxy surveys, which find that the disk 
galaxy population appears to have undergone little 
growth since z ~ 1 (e.g., iMelbourne et all 120071 : 



Sargent et all 120071: IPapovich et all 12003 jBarden et all 



2003 iRavindranath et al J l2004lLiUveta!]ll99fij ). While 
M81 is only one large disk galaxy and our field is only 
a small portion of it, the similarity to Milky Way, M31, 
and even lower mass M33's stellar populations supports 
the same scenario. Measurements of the stellar popu- 
lations of local disks therefore strengthen the results of 
some galaxy surveys by attacking the problem with a 
completely independent technique and finding a similar 
result. 

Finally, the spatial distribution of main-sequence stars 
in the field show that the field partially covers the out- 
skirts of a spiral arm as well as an interarm region. De- 
tailed analysis of the resulting CMDs for the contrasting 
regions shows that the difference between the popula- 
tions is due to the fraction of stars with ages younger 
than ~0.5 Gyr ago and results from only a small per- 
centage (~2%) of the stars, similar to the arm and inter- 
arm p opulations seen in M31 disk populations (jWilliamsl 
2002). In addition, stars younger than ~100 Myr appear 
to be confined to the arm region, suggesting the structure 
survives for at least 100 Myr. 

Support for this work was provided by NASA through 
grant GO-10915 from the Space Telescope Science Insti- 
tute, which is operated by the Association of Universities 
for Research in Astronomy, Incorporated, under NASA 
contract NAS5-26555. 
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TABLE 1 

Distances and Extinction Values from MATCH 



Subficld (m - M) A v 



Full 27.89±0.05 0.38±0.06 

Arm 27.88±0.08 0.45±0.06 

Interarm 27.92±0.06 0.33±0.06 

Crowd 27.88±0.06 0.32±0.06 



TABLE 2 

Cumulative Fraction of Stars Formed at Each 
Age 



Redshift Age (Gyr) Fraction of Stars Formed 

0.001 
0.002 
0.004 
0.006 
0.010 
0.015 
0.024 
0.049 
0.216 
1.080 



0.01 


i.oooIq 


0.03 


0.999tg 


0.05 


0.999t° 


0.08 


0.999^0 


0.13 


0.999lo 


0.20 


0.997to 


0.32 


0.994lg 


0.63 


0.979±" 


2.51 


0.848lo 


7.94 


0.622^^ 


14.1 


o.oootS 
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Fig. 1. — The locations of our M81-DEEP field and our defined subregions of the field shown on a DSS image. The arrow marks the 
apparent spur of the Northern spiral arm according to the distribution of main-sequence stars in our fiel d. Inset on the lower- right corner 
is the equivalent location of our M81-DEEP field shown with a white box on a star count map of M31 (Ferguson ct al. 2002). The inner 
and outer edges of the field are located at the same scale lengths as in M81. 
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Fig. 2. — Histogram: The distribution of the combined sharpness values for objects in our initial photometry catalog. Vertical Line: The 
sharpness cut we applied to our final photometry catalog. 
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Fig . 3. — Left: The F606W, F814W CMD for our entire ACS field. Lines show a small subset of the isochrones used to fit the 
data (Marigo ct al. 2008, shifted assuming Ay = 0.3 and (m — M)q = 27.9). Isochrones shown are (from blue to red): [M/H]=-0.4, 
log(age)=7.3; [M/H]=-0.4, log(age)=7.6; [M/H]=-0.4, log(age)=8.0; [M/H]=-0.4, log(age)=8.3; [M/H]=-0.4, log(age)=8.6; [M/H]=-1.3, 
log(age)=10; [M/H]=-0.7, log(age)=10; [M/H]=-0.4, log(age)=10; [M/H]=-0.2, log(age)=10; [M/H]=0, log(age)=10. Middle: The CMD 
of the arm region shown in Figure [4] Right: The CMD of the interarm region shown in Figure [4] In areas where the points would saturate 
the plot, we provide contours following the density of points in that part of the CMD. 




Fig. 4. — The borders of the subregions used in our SFH analysis. Blue crosses mark the locations of main sequence stars in the field. 
Main sequence stars were chosen using a hand-drawn polygon that followed the edges of the blue plume of stars in the CMD. The color 
and magnitude limits of the polygon were approximately 24<F814W<28 and -0.2<F606W-F814W<0.3. 
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Fig. 5. — Top: The mean residual magnitude and root-mcan-square error of the artificial star tests in the F606W filter measurements 
are shown as a function of input star magnitude. Bottom: Same as Top, but for the F8f4W filter measurements. 
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Fig. 6. — The completeness measured from our artificial star tests is shown in grayscale as a function of color and magnitude. The scale 
is linear, with 100% completeness shown as black and 0% completeness shown as white. White diamonds mark the completeness limit used 
to fit the SFH of the field. 
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Fig. 7. — Left: Gaussian plus line fits to our measured magnitude histogram to determine the magnitude of the red clump (drawn in 
red) and AGB bump (drawn in blue). Right: Gaussian fits to the color histograms measured at the best-fit magnitude of the red clump 
(red) and AGB bump (blue). 
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Fig. 8. — Left: Small portion of our full-field CMD centered on the red clump region. The black circle marks the best-fitting center for 
the feature. The white circle marks the measured height and width of the feature _„„„„. „ ,„,i = 0.154; cr m „„. =0.191). Boxes 
mark the locations of the fits to the same feature in model CMDs convolved with our photometric errors and completeness. Redder colors 
denote higher metallicities (the metallicity range is -1.0<[M/H]<0.2); larger boxes denote older ages (2 Gyr<age<13 Gyr). Large error 
bars in the upper right show the 1-sigma width of the feature in a single age and metallicity CMD ( CT (m F606w -m F814w ) = 0.129; <r mFgl4w 
= 0.146). Arrows show our 1-sigma error ranges for distance and reddening. (Right: Same as Left, but for a small portion of our full-field 
CMD centered on the AGB bump region. The measured height and width of the feature are CT (m F60 6w- m F8i4iv) = 0-148 and fm Fsl4w = 
0.123. Values for the 1-sigma width of the feature in the a single age and metallicity CMD are CT (m F606w - — m F814w ) = 0.070 and o" mF814vv 
= 0.073) Note the fitted AGB bump center is slightly brighter than the mode value due the the slope of the luminosity function skewing 
the mode to fainter magnitudes. Predicted magnitudes assume Ay = 0.3 and (m — M)q = 28.0. 
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Fig. 9. — The age-metallicity plane. Contours denote different Nj\qb^ : Nrc ratios, increasing from 0.026 to 0.051 from left to right 
and labeled. The shaded area denotes where model stellar populations have a ratio N^GBb '■ Nrc consistent with that measured from our 
data. 
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Fig. 10. — Our best full CMD fit to the data from the entire field. Upper-left: The observed CMD. Upper-right: The best-fitting 
model CMD from MATCH. Lower-left: The residual CMD. Redder colors denote an overproduction of model stars. Bluer colors denote an 
underproduction of model stars. Lower-right: The deviations shown in lower-left normalized by the Poisson error in each CMD bin. This 
plot shows the significance of the residuals in lower-left. Only the red clump and AGB bump show significant residuals. 
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Fig. 11. — Histograms of the luminosity functions of our best full CMD fit to the data from the entire field. Black: The observed 
luminosity function of the entire field. Gray: The luminosity function of the best-fit model CMD. 
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Fig. 12. — The SFH of the entire ACS field as determined by the MATCH package. Top: The solid histogram marks the star formation 
rate (normalized by sky area) as a function of time for the past 14 Gyr. The dashed line marks the best-fitting constant star formation rate 
model. Middle: The mean mctallicity and metallicity range of the population as a function of time. Heavy error bars mark the measured 
mctallicity range, and lighter error bars mark how that range can slide because of errors in the mean metallicity. Bottom: Same as top, 
but showing only the results for the past 1.3 Gyr. 
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Fig. 13 . — The SFH of the arm subregion as determined by the MATCH package. Lines, error bars, and panels are the same as in 
Figure [12] 



22 



Redshift (z) 

52 1 0.5 0.2 0.1 0.05 0.02 0.01 0.005 0.002 0.001 




10.00 



1.00 0.10 

Age (Gyr) 



0.01 



Redshift (z) 

52 1 0.5 0.2 0.1 0.05 0.02 0.01 0.005 0.002 0.001 




10.00 



1.00 0.10 

Age (Gyr) 

Redshift (z) 
0.05 



0.01 



0.02 0.01 




0.8 0.6 0.4 

Age (Gyr) 



0.0 



Fig. 14. — The SFH of the interarm subregion as determined by the MATCH package. Lines, error bars, and panels are the same as in 
Figure [12] 
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Fig. 15. — The SF H of the crowded southern subregion as determined by the MATCH package. Lines, error bars, and panels are the 
same as in Figure [12] 
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Fig. 16. — The SFH of all subregions are overplotted. Rates are not normalized by area in order to show the contribution of each 
subregion to the star formation of the total field. Black: The SFH of the full field as determined by the MATCH package. Blue: The SFH 
of the arm subregion as determined by the MATCH package. Green: The SFH of the crowded, southern subregion as determined by the 
MATCH package. Red: The SFH of the interarm subregion as determined by the MATCH package. 
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Fig. 17. — Black: The normalized cumulative star formation of the full field as determined by the MATCH package. 
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Fig. 18. — Met allicity range of the dominan t stellar populations in deep resolved photometry for M 31 l^rowri etaDHSSD, M33 
flBarker et al.ll2007ft , M81 (this work), NGC 5128 JReikuba et al.|[2005f) . and NGC 3377 JHarris et ai.ll2007T) are plotted along with that of 
the Milky Way thick disk llAllende Prieto et al.||2 006l) a gainst several other properties. Top: The metallicities as a function of the absolute 
K- band magnitude of the galaxy (Skrutskic ct al. 2003). The solid and dashed lines show t he lum inosity-metallicity relations determined 
by Trcmonti ct al. ( 2004j, gas phase metallicity of galaxy central regions) and Mouhcinc (20061, stellar red peak metallicities of galaxy 
"halos"), respectively. These relations were converted from B-band and V-band to K-band using the Tully-Fisher calibrations of IVerheiienl 
(2001) and Sakai ct al. (20Q0j). The Milky Way luminosity was calculated by applying V ro t= 220 km s 1 to the Tully-Fisher calibration of 
IVerheii cn (2001). Middle: The populations' metallicities as a function of the radii at which they were sampled (normalized to the half-light 
radius of the galaxy). Bottom: The populations' metallicities as a function of their ages. 



